Degrees of Freedom in Massive Gravity 



D. Comelli'^, M. Crisostomi'''^ F. Ncsti'' and L. Pilo''^'= 
''INFN - Sezione di Ferrara, 1-35131 Ferrara, Italy 
^ Dipartimento di Fisica, Universitd di L'Aquila, 1-67010 L'Aquila 
'^INFN, Laboratori Nazionali del Gran Sasso, 1-67010 Assergi, Italy 
(Dated: February 20, 2013) 

We study in a systematic way a generic non derivative (massive) deformation of general relativity 
using the Hamiltonian formalism. The number of propagating degrees of freedom is analyzed in 
a non-perturbative and background independent way. We show that the condition of having only 
five propagating degrees of freedom can be cast in a set of differential equations for the deforming 
potential. Though the conditions are rather restrictive, many solutions can be found. 
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Whether general relativity (GR) is an isolated theory 
is an interesting question both from the theoretical and 
phcnomcnological side. In particular, attempts to modify 
GR at large distances recently have received a renewed 
interest. In this paper we focus on massive deforma- 
tion of GR. A possible definition of Lorentz-invariant (LI) 
massive gravity is any deformation of GR such that the 
metric perturbation /i^i, around flat Minkowski space 77^1, 
behaves as a massive spin 2 particle with 5 polarizations 
(DoF) in contrast with GR where only 2 DoF propagate. 

Since the work of Fierz and Pauli (FP) [l[ it was re- 
alized that for a generic massive deformation, even at 
linearized level 6 DoF are present. Though a tuning al- 
lows to get rid of the sixth mode, it reappears at the 
non-linear level or around non flat backgrounds [Ij. In 
addition the sixth mode is typically a ghost. Recently, it 
was shown that there exists a LI non linear completion 
of the FP theory that is free of ghost up to the fourth 
order Q avoiding the presence of the sixth mode. Then 
proof of the propagation of only five DoF and the absence 
of instabilities was extended to the non perturbative level 
in 0, IBl ■ Cosmology of such a model was studied in @ . 
On the other hand, one can give up Lorentz invariance in 
the gravitational sector, and a class of theories of massive 
gravity that avoid the vDVZ discontinuity [?[ is available 
at quadratic level (8l-[ll|. 

The main goal of this paper is to study the number 
of propagating modes in generic massive deformations of 
GR, non perturbatively and in a background indepen- 
dent way. We stress that our analysis applies to both the 
Lorentz invariant and Lorentz breaking case. For a dif- 
ferent approach that focus only on the Lorentz invariant 
potential in 0] see [T2| . 

The natural tool for such a task is the Hamiltonian for- 
mulation in the Arnowitt-Descr-Misncr (ADM) [l^ for- 
malism. The 10 DoF of the metric are split in the lapse 
N, the shift iV' and the spatial 3-metric jij (i,j=l,2,3) 



(1) 



Any massive deformation of gravity can be obtained by 
adding to the Einstein-Hilbert action an arbitrary non- 
derivative function V{N^ N'\^ij) of these variables. Thus 



the Hamiltonian reads 



H = M^i / d-'x [N"^ Ha + N ^V] 



(2) 



where we introduced the notation N"^ ^ Na = {N, N') 
and Ha = {H, Hi), with A = 0,1,2,3. The derivative 
dF/dN^ of a function F with respect to will be 
denoted by F,a, and we define V = m'^N-f^l'^V, with 
7 = det(7y ). As matter of fact. Ha are functions of the 
spatial metric 7^- and its conjugate momenta H'-' only. 

When m — 0, the action ([2]) reduces to the familiar 
ADM GR Hamiltonian derived from the Einstein-Hilbert 
action. In general, this action can describe Lorentz in- 
variant (LI) or Lorentz breaking (LB) models [1, 

Though V in ([2]) breaks diffeomorphisms, by adding 
a suitable set of Stiickelberg scalar fields gauge invari- 
ance can be restored IJ] and ([2]) can be considered as a 
generally covariant theory written in a unitary gauge. 

Peculiar of GR and of its deformed version ([2|) , is that 
are non dynamical and that their conjugate mo- 
menta Ha are vanishing on the physical phase space, 
which are called primary constraints. The canonical 
analysis of constraints is thus appropriate to assess non- 
perturbatively the number of degrees of freedom DoF. 

Canonical analysis. The total Hamiltonian is 



H 



d^x V'Ua ■ 



(3) 



where a set of Lagrange multipliers A"^ are introduced 
to enforce the primary constraints H^ w 00 The time 
evolution of any function of 7^ , and their momenta 
is given by the Poisson bracket with Hq^: 

df{t,x) 



dt 



-{f{t,x),HT{t)]^{f{t,x),H{t)} 



d'y\^{t,y){F{t,x),nA{t,v)}. 



(4) 



The Poisson bracket between two functionals F and G of 
canonical variables qa , p° is defined as 

6F 5G SF SG 



{F, G)^j d^y 



6qa{t, y) Sp'^it, y) Spa{t, y) 5q''{t, y) 



^ The customary notation 11^ K, means that Ha arc only weakly 
zero, e.g. zero on the constraints surface only. 
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In the following the dependence on space and time of the 
various fields will be understood. 

Even before studying the equations of motion of the 
dynamical variables hij, it is crucial to enforce that the 
four primary constraints are conserved in time. This 
leads to secondary constraints Sa- 



Sa = {TIa,Ht} = -{Ha + V.a) ~ . 



(5) 



Basically, these are the equations of motion of N^. They 
do not determine any of the Lagrange multipliers and are 
not automatically satisfied. Therefore, Sa ~ reduce the 
dimension of the physical phase space by 4. Also Sa must 
be conserved, leading to the new conditions 

Ta = {Sa, Ht} = {Sa, H} - Vab 



0. 



(6) 



Vai 



V 



ab 



dN^dNB • 

Now, if the Hessian \Vab\ is non-degenerate, it can be 
inverted and the four eqs. ([6]) can be solved for all the 
A^. At this point all constraints are consistent with time 
evolution and the procedure stops. 

Summarizing, in case det \Vab\ 7^ we have 4 (11^) + 
4(5^) = 8 constraints, for a total of (20 - 8)/2 = 6 
DoF. Technically, these 8 constraints are Second Class, 
being Rank|{nyi, Sb}\ = R.ank|VA_B| =4. As a result, 
the whole gauge invariances of GR is broken. 

In the LI case, expanding around a Minkowski back- 
ground, Lorentz invariance tells us that the propagating 
DoF must be grouped in a massive spin two (5 DoF) rep- 
resentation plus a scalar (1 DoF). This is the Boulware- 
Deser result, valid for a generic potential. The extra 
scalar, the so called Boulware-Deser sixth mode 0], turns 
out to be a ghost rendering a generic nonlinear comple- 
tion unviable. 

The first important result is thus that a necessary con- 
dition to have less than six propagating DoF, is that 

def 

r = RanklVAsI < 4. 

The case r ~ 3. In this important case the matrix Vab 
has a null eigenvector and three other eigenvectors 
— 1' 2, 3) with eigenvalues k„: 



Vab x' 



0. 



Vab = K„ Ei 



(7) 



For instance if dei{V^j) ^ 0, then x'^ = (1, -V,:]^ Vqj). 
The Lagrange multipliers can also be split along and 
its orthogonal complement 



\^ = zx^ + 



3 

E 



ZX^ + A" 



(8) 



The Hessian is now non-degenerate only in the sub- 
space orthogonal to x^i and from ([SJ one can deter- 
mine only three out of four Lagrange multipliers, d„ = 
E^{Sa, H}/Kn- The projection along x^ leaves z unde- 
termined and ^ gives a new non-trivial (tertiary) con- 
straint 



r^x 



^{Sa,H} 



0. 



(9) 



Here and in the following we suppose that the constraints 
found are non trivial. At this point, either the conserva- 
tion of T allows to determine the remaining multiplier z 
or it may generate a further constraint. One has 

Q^{T,Ht} = {T,H} + {T,X^ -Ua}, (10) 

\^-nA = Jd^yX^{y)UA{y). 

Clearly, the first piece of (|10p does not contain z. Using 
([7]) and ^ after a tedious but straightforward computa- 
tion we find for the second piece 



{T, A^ • H^} « x^{X^ ■ Sa,Sb} - X^{A^Vba, H} 

, ax^ 



A^PVsD-e-z, 



(11) 



e{x, y) = x^(x) {Sa{x), Sum x'^iy) ■ (12) 

It turns out that Q[x,y) = A^{x^y)di5^'^\x — y) with 
A(x, y) = A(y, x); it is worth to point out that in a finite 
dimensional context would be absent. Thus, 

e • 2 = jd'y e(x, y)z{y) = -^^5, [z{xfA\x, x)] . 

(13) 

As a result, Q is free from z \i A^{x, x) — 0, which consists 
in the following new condition 



X 



°^V. + 2xV^ 



0, 



V = -f-^/'^V. (14) 



If this condition is satisfied, Q in pOj) is a quaternary 
constraint, which together with T confirms the elimina- 
tion of one DoF. If condition ([T4)) is not satisfied, ([T0| is 
not a constraint but can be used to determine the last 
Lagrange multiplier z. This leaves us with an odd num- 
ber of second class constraints and an odd dimensional 
physical phase space, in particular 54-1/2 DoF. A similar 
situation is encountered in other deformations of GR, see 
for instance . It not clear if such a peculiar occurrence 
is physically acceptable or not. For sure, condition 
is rather restrictive. 

When p4|) is satisfied by V, the time evolution of the 
new constraint Q is non trivial and can finally determine 
the multiplier 



{x^Q.a) [A^Q,B + {Q,ff}] 



(15) 



In conclusion, for r = 3, there are 4(H^) + A{Sa) + 
1 (T) + 1 (Q) = 10 constraints, giving (20 - 10)/2 = 5 
DoF, that is a good candidate for a theory of a massive 
spin-2. The required conditions arc that holds and 

Rank|VAB| = 3 and x^ Q,a 7^ ^ DoF = 5. (16) 

To our knowledge, the observation that a singular Hes- 
sian is a necessary condition to avoid the Boulware-Deser 
argument which leads to six propagating DoF was made 
for the first time in general terms in Q . 
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We stress that the condition p4|) does not appear at 
the hnearized level around flat space; indeed, the depen- 
dence on z starts at the quadratic level thus it manifests 
either at the non-linear level or in a non trivial back- 
ground. As a result, a random extension of the non-linear 
level of the quadratic action of Q will not satisfy (fT4l). 

General case. Let us briefly discuss the generic case 
with the Hessian being of any Rank|Vyi_B| = r between 
zero and four. Now Vab has 4 — r null eigenvectors Xa 



Vab xf = 0, 



1, 



(17) 



and as before, we can split the Lagrange multipliers as 



4— r 

E 



(18) 



The projections of © along determine the r La- 
grange multipliers dn = E^{Sa, H} / Kn, while the pro- 
jections along the Xa S^^^ 4 — r tertiary constraints 



The conservation of these contraints leads to 



(19) 



Qa = {%,H} + Y,{Ta:nA}-dnE'„ 



4— r 
13 = 1 



e^p ■ zp (20) 



which are 4 — r relations linear in the remaining 4 — r 
Lagrange multipliers Zq. The DoF thus depend on how 
many multipliers are determined in these conditions. The 
matrix Oap is defined as 

9^0 = xiit,x){SAit,x),SB{t,y)}xpit,y). (21) 

Notice that due to the non-trivial dependence on x, y, 
6ap is not necessarily antisymmetric and its rank s is not 
always even. If s < 4 — r, some Za arc undetermined and 
we have 4 — r — s new quartic constraints which reduce 
the DoF. 

At the next (fifth) step, if the conservation of quar- 
tic constraints determines all the needed multipliers, the 
procedure stops and we have only the 16 — 2r — s con- 
straints 4 (n^) + 4 (5a) + (4 - r) (To) + (4 - r - s) (Q„). 
This implies that the maximal number of states is 



DoF < 4 + r/2. 



(22) 



This is the maximal number because if e.g. some of the 
Za are not determined further steps are necessary, pos- 
sibly reducing the number of DoF. In general also first 
class constraint may be present corresponding to residual 
gauge invariances, further reducing the number of DoF. 

Summarizing, massive deformations with 5 DoF exist 
only in two cases: r = 3, s = and r = 2, s = 2. In 
the case r = 3, the condition s = is equivalent to the 
differential equation p4[) for the deforming potential. In 
the case r = 2 and s = 2, (j20p determine the remaining 



two LMs and all 10 constraints are consistent with the 
time evolution. Notice that it is not possible to build a 
potential with 5 = 2 without breaking spatial rotations. 
Thus, we have shown that r = Rank(VAB) = 3, together 
with (fT4)) are necessary and sufficient conditions for a 
rotational invariant massive gravity theory with 5 DoF 
(again we suppose that none of the found constraints are 
trivial, i.e. ~ 0). 

Solutions. Let us then focus on how the r = 3 condition 
can be realized. If spatial rotations are preserved by the 
potential, the only option is 



det{VAB)=0, dct(V.y)^0 



(23) 



or equivalently Voo — Voi(Vij)^^ Vjo = 0. Physically this 
means that once the equation for the shifts N"^ are solved 
and inserted back in the potential, V will be linear in N. 
When and are satisfied, V is a candidate for 
a theory of massive gravity with 5 propagating DoF. It 
is well known that pcrturbatively the number of propa- 
gating DoF depends strongly of the background [131 ■ On 
the contrary, our analysis is fully non pcrturbative, back- 
ground independent and it is not sensitive on the way the 
deformation is built. For instance, it equally applies to 
the case where V is built from the invariants of the ma- 
trix f, where / is an arbitrary fixed metric or is 
built out of set of four scalar fields [1, Q . 

Let us look more closely to ((23)) . As a general remark, 
it is a homogeneous Monge- Ampere equations [Tcj ; a par- 
ticular solution is known in a closed form and, remark- 
ably, also the general solution can be given in paramet- 
ric form. Thus, there exist a large class of potentials, 
taken as a function of lapse, shift and spatial metric 
V{N, N\ 7ij ) that satisfies Typically, one first finds 
a family of solutions of ([25]) and then tries to enforce (|14p . 

In fact, it is not hard to find particular solutions of 
(1231), though we stress that most of the solutions found 
will in general be LB. For instance, as a simplification 
one can take V to be function of just two rotationally 
invariant combinations of lapse and shift, x = N, y = 
jijN^N^, i.e. V = V{x,y). In this case ([^5)) reduces to 



2y 



VyV, 



(24) 



for which an infinite class of solutions can be found by 
separation of variables. An example is a class of solutions 
that we will compare later with the LI ones 



V = /3i {x + - iy'/-' + Ps) 



1/2 



+ f^iX, 



(25) 



where /?«=!,. ..,4 are arbitrary scalar functions of 7ij. Thus 
it is relatively simple to construct solutions of . Much 
more restrictive is condition (|14[) . In two dimensions, 
the class of potentials ([25)1 satisfies also (ITil) if f3i = ci, 
(32 = C2"/^^^, Pi = C47"'^/^ and fS^ is an arbitrary function 
of the spatial metric 711, with ci, C2, C4 constants. 
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More in general, it is known that in four dimensions 
the general solution of the Monge- Ampere equation can 
be written implicitly in terms of two functions [igI ]. 
When (fn)) is also considered, one is left with infinitely 
many solutions in terms of a single function. Of course 
most of the solutions will be Lorentz breaking. An ex- 
ample is the following two-parameter solution 



where g^-' are the spatial components of the inverse 4d 
metric. One can check that both ([25]) and ([H)) are satis- 
fied. 

Generic LB solutions with five propagating DoF are at- 
tractive in their own respect. Indeed, expanding such a 
LB potential around flat space at quadratic level, one can 
easily get a theory featuring a finite range (Yukawa-like) 
gravitational potential @, The breaking of Lorentz in- 
variance, taking place only in the gravitational sector, is 
not phenomenologically excluded (and on the contrary is 
potentially observable in the forthcoming gravitational 
waves experiments). More interestingly, the resulting 
theory can be free from the vDVZ discontinuity and re- 
main perturbative at the solar system scale. Note that 
on the contrary to avoid the vDVZ discontinuity an y L I 
theory needs to rely on the Vainshtein mechanism |18j . 
which is not yet completely understood. Even taking 
it for granted, the price to maintain LI is high: grav- 
ity in the solar system is non-perturbative, making the 
theory hardly predictive. Thus Lorentz breaking seems 
to offer a healthier alternative. After all in the presence 
of sources Lorentz symmetry in the gravitational sector 
plays little role if any. For instance, to check LI in the 
gravitational sector one should measure the velocity of a 
graviton in (almost) flat space and compare it with the 
speed of light. We leave for a future work the complete 
study of the solutions of the above differential equations 
leading to 5 DoF and of their phenomenology. 

LI solutions. Let us now study in detail the problem of 
finding a ghost-free and Lorentz invariant potential. The 
requirement of Lorentz invariance reduces considerably 
the possible solutions of (|23p . We require that the poten- 
tial is built out of the invariants of X = g~^'rj, made out 
from the spacetime metric and a fixed Minkowski metric 
77, which automatically guaranties that V is Lorentz in- 
variant. To start, let us consider the simplest case of two 
dimensions. The potential V is a function of the spatial 
metric 7 = 711, the lapse N and shift N'^, and eq. 
reads 



V 



0. 



(27) 



To single out the LI solutions, one observes that V = 
V{X) is only function of the eigenvalues Ai_2 of the matrix 
X. Therefore, after expressing TV, in terms of Ai, A2, 
7, (|27p must hold for any 7. The resulting equation is 
cubic in it, so that ([?7|) splits into two branches of three 



differential equations 
3 



i;(2,0) 



V 



2Ai ' 2A2 



2AiA2(A;/'±A^/') 



(28) 



(26) Both can be solved in terms of simple functions: 



>;1a/AiA2 + Q!2 (\/A7± VA^) + Ct3 

V A1A2 



(29) 



with Q!i,2.3 arc integration constants. These solutions also 
satisfy . Indeed, they are special cases of (^5]) , as can 
be seen by writing them as functions of x, y, 7: 



1/2 



Being LI, they can also be written as functions of X 
Tr(Xi/2) a. 



Vi ~ ai + a2 ■ 



Vii ^ ai + a2 



VdctX VdctX 



Ti{X) - 2VdctX 



"3 



VdBtX 



VdctX 



(30) 



(31) 



The first is precisely the two-dimensional version of the 
ghost-free potential found in Q. The solution corre- 
sponding to the — sign on the other hand is a new one. 
We stress that the two family of solutions do not overlap. 
It can be shown that the theory with Vii can not admit 
a Minkowski background, being non-analytic there. On 
the other hand Vi does. 

The same strategy can be extended to d = 3 and 4, 
though the computation becomes much more involved. 
For d ~ 3 from (|23|) we obtain a system of 9 partial dif- 
ferential equations which admit four families of different 
solutions, generalizing ([29|) to the case of three eigenval- 
ues, with possible switching of signs. In four dimensions 
there are eight solutions, which satisfy a system of 125 
partial differential equations. The first solution coincides 
with the one found in Q, generalizing ([5(7)) to d = 4. The 
other solutions are again obtained by permuting the signs 
of the eigenvalue square roots, and generalize ([51]) . Also 
in the four dimensional case only V/ admits a Minkowski 
background. This is not very surprising, as the solution 
in [3| was obtained after a resummation of the pertur- 
bative theory around Minkowski. Due to the compli- 
cated coupled system of partially differential equations, 
in d = 3 and d = 4 we cannot prove that the solutions 
found in the LI case are unique, though the 2d analysis 
seems to suggest so. In d = 3,4 one can show that the 
generalization of (|30p satisfy also p3 Comparing the LI 
with the LB cases it is clear that Lorentz invariance is 
a severe constraint that cuts off most of the solutions of 
even before dealing with ([H]). 
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Discussion and conclusions. We studied the Hamil- 
toniaii structure of a general massive deformation of 
GR; which leads to non-perturbative and background- 
independent results. In particular, the analysis applies 
to both the LI and LB cases. The condition for having 5 
DoF amounts to a set of differential equations for the de- 
forming potential: (j23|) . which is of Monge- Ampere form, 
and (113]). The picture that emerges is that Lorentz in- 
variance considerably restricts the potentials with 5 DoF: 
in the lower dimensional case ci = 2, a unique solution ad- 
mits a Minkowski background. In the d = 3, 4 we checked 
that the ghost free potential proposed in [H satisfies ([25)1 
and (fT4|) . In general there are new LI potentials, which 



however do not admit Minkowski as a background. Phe- 
nomenologically, LI theories are hard to handle due to 
their non-perturbative nature inside the Vainshtein ra- 
dius, that has to be larger than solar system to avoid the 
clash with the bending of light experiments. As a result, 
the theory is highly nonlinear and even predicting the 
motion of Mercury becomes non-trivial. On the other 
hand, in our investigation we have clearly shown that 
there are many Lorentz breaking solutions of (|23|) and 
(fM)) which are candidates to be a consistent and weakly 
coupled phase of massive gravity. The analysis of the so- 
lutions of and ([T^ and their phenomenology will be 
given in a future work. 
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